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We discuss characterization of the tuned-harmonic mode of operation in 
IMPATT oscillators, and introduce an equivalent circuit which incorpo- 
rates the large-signal, " 'single-frequency" oscillator admittances at the 
fundamental and second-harmonic frequencies. Complete characterization 
of this mode is equivalent to specifying the behavior of each of the four 
elements of the equivalent circuit as functions of the oscillation state vari- 
ables: fundamental voltage and frequency, second-harmonic voltage and 
relative phase. Using the approximate large-signal analysis of Blue, 1 the 
values of the equivalent circuit elements are presented, as an example, for 
a 6-GHz IMPATT diode under a variety of oscillation conditions. This 
equivalent circuit is used to clarify the role played by the fundamental and 
second-harmonic, single-frequency oscillator admittances in the tuned- 
harmonic mode. 

Using an approximation to the equivalent circuit, we investigate the 
criteria for stable oscillation of the tuned-harmonic mode. It is found that 
the stability criteria are in general quite restrictive. For the same 6-GHz 
germanium diode, the range of stable phase is investigated, as a function 
of the RF parameters, for certain special cases. It is found to be possible 
to satisfy the stability criteria for the phase which gives an optimum en- 
hancement of the fundamental power output if certain conditions on the 
external RF circuit are satisfied. 

I. INTRODUCTION 

It was found by Swan 2 that the introduction of a trapped resonance 
at the second harmonic of the oscillation frequency in a 6-GHz Ge 
IMPATT diode oscillator provided dramatic increases in the output 
power and efficiency, as compared with the results obtained with the 
ordinary single quarter-wave transformer coaxial circuit. Since that 
time several authors 1 ,3_s have reported both theoretical and experi- 
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mental examinations of the effect. It appears that the addition of a 
properly phased second-harmonic voltage improves the phasing of the 
RF current relative to the fundamental voltage so as to increase the 
negative conductance and (at least at lower frequencies) give an in- 
crease in the power output at the fundamental frequency. The circuit 
conditions required for the observation of this effect have been incom- 
pletely understood. 

The purpose of this paper is to present the results of an analytical 
study of the interaction of an IMP ATT diode with a circuit having 
resonances at two harmonically related frequencies. The analysis is 
begun by the introduction of an equivalent circuit for the diode by 
which these two-frequency oscillators may be characterized. A stability 
theory is then developed along the lines taken by Kurokawa which 
examines whether a particular circuit, even though matching the 
impedances required by the diode at both frequencies, will or will not 
provide a stable oscillation. 9,10 The stability theory is examined in 
some generality, and three special cases are studied for which tractable 
analytical results can be obtained. It is found that in the case of zero 
fundamental or second-harmonic voltage, the theory reduces to the 
single-frequency stability criteria derived by Kurokawa. In more 
general cases, the theory indicates that by designing (or adjusting) 
the circuit carefully one can obtain stable operation at phase angles 
which enhance the fundamental power. However, the theory also 
indicates that stable operation may be impossible if the circuit-diode 
interaction is not just right, even though the diode and circuit are 
matched to each other at the two frequencies. 

In a final section, a numerical example is given in which the theory 
is applied to a model of a 6-GHz germanium IMPATT diode, using 
the approximate large-signal analysis of Blue. 

II. TWO-FREQUENCY CHARACTERIZATION 

The IMPATT oscillator is truly a single-frequency oscillator only at 
very small ac voltages and currents. At larger signal levels the non- 
linearity is very strong, and therefore there should be strong inter- 
actions between harmonically related signals. However, by operating 
the diode in a well-designed single-frequency circuit, the power output 
may be limited to a single frequency. This may be done, for example, 
by presenting short-circuit, open-circuit, or reactive loads at the har- 
monic frequencies. In the case of short circuited harmonics, the har- 
monic voltage amplitudes V k are zero, and only the fundamental 
voltage Vi is nonzero. It is then common practice to calculate a large- 
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signal diode admittance as a function of V 1 and to use this admittance 
to describe device behavior. On the other hand, for the case of open- 
circuited harmonics, the harmonic currents, I k are zero, and only the 
fundamental current 7 : is nonzero. It is then preferable to characterize 
the diode by a large-signal impedance which is a function of the RF 
current amplitude I x . Both of these conditions constitute tunings at 
the harmonic frequencies, albeit ones that are particularly useful and 
simple to express analytically. 

To consider other, more general, loading conditions at the harmonic 
frequencies, one must introduce two more variables (amplitude and 
phase) for each additional frequency for which the amplitude is nonzero. 
One of the most important points is that the input admittance (for 
example) at the fundamental frequency is no longer a unique function 
of Vi and the frequency /; but instead denning the state of oscillation 
requires a vector whose components are V x , • • • , V N , f, cp 2 , ■ ■ ■ , <pn 
where N is the maximum harmonic number of interest and <p k is the 
phase of the fcth harmonic voltage relative to the fundamental. This 
vector does uniquely describe the state of oscillation, and for every 
such vector, there exists a set of complex admittances iji • • • y N which 
are uniquely determined. If this is not so, it simpty means we have 
inadequately described the system and must include more component 
signals, either harmonics or subharmonics. 

We shall limit the discussion to include only two harmonically related 
frequencies and consider that V k = for k > 2. This also means that 
we will only discuss the admittance characterization and not the im- 
pedance characterization. 

A convenient way of utilizing the information already known about 
the large-signal single-frequency admittance of the diode is to separate 
the input admittances at the two frequencies as shown in Fig. 1. This 
equivalent circuit shows a fundamental port and a second-harmonic 
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Fig. 1 — Equivalent circuit of the IMPATT diode which includes nonzero voltages 
at two harmonically related frequencies. Port 1 is the fundamental port and port 2 
is the second-harmonic port; y n and ?/ 22 are the large-signal single-frequency diode 
admittances at the fundamental and second-harmonic frequencies, respectively. 
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port. The admittances y n (Vi) and y 22 (V 2 ) are the large-signal single- 
frequency admittances that would be measured at the fundamental if 
there were no harmonic (or subharmonic) voltages present. That is, 
they are just the ordinary large-signal admittances y(V) at the fre- 
quencies / and 2/. 

The admittances y 12 0\ , Y* , /, <p 2 ) and y n (Vi , V 2 , f, <p 2 ) account 
for the conversion of current between the two frequencies and it is the 
study of their effects that is the main subject of this paper. The phase 
<p 2 is defined by the assumed voltage waveforms 

"i(0 — Vi cos u) n l 
and 

v 2 (t) = V 2 cos (2u) n t + tp 2 ). 
The input admittances are 

„ , Va exp Q 2 ) m 

Y inX = y u + ijy. y (1) 

and 

>••,„, = fe + V« 7 , a 5 eW ® 

at the fundamental and second-harmonic frequencies respectively. 
Since y u and y 22 are independent of the phase <p 2 by definition, equations 
(1) and (2) show that the input admittance loci for fixed V\ and V 2 will 
be counter rotating closed curves as a function of <p 2 . These curves will 
enclose the admittance points y n and y 22 separately providing that 
y l2 and y 2 \ are not strong functions of ip 2 . If, for example, y l2 and y 2l 
are independent of <p 2 , Y inl and Y in2 will be circles centered about 
y n and y 22 respectively, the radii of which depend upon the ratio 
V 2 /V x . They generally turn out to be somewhat elliptical in shape 8 
although, in many cases, of very low eccentricity. 

Figure 2 is the calculated 1 large-signal, single-frequency, admittance 
plane plot for a 6-GHz germanium diode, from which y n and y 22 may 
be obtained directly. Figures 3 and 4 show Y inl and F, n2 for various 
fundamental frequencies when the voltages are held constant, demon- 
strating the elliptical and circular behavior noted above. Note that in 
Fig. 4 the second-harmonic input admittance has a positive real part 
for some ranges of the phase <p 2 . To operate at such phase angles and 
RF voltages, the external circuit must supply power to the diode at 
the second-harmonic frequency, and thus these conditions are un- 
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Fig. 2 — The calculated large-signal single-frequency admittance of a 6-GHz 
germanium IMPATT diode at a bias current density J = 340 A/cm 2 . 



realizable when operating into a passive circuit. The diameter of these 
admittance contours is inversely proportional to the second-harmonic 
voltage amplitude V 2 , however, so that at higher values of V 2 , the 
entire contour may lie in the left-half plane. 

The rather simple structure of the Y inl and Y in2 loci of Figs. 3 and 4 
suggests that y l2 and y 2i might be rather insensitive functions of <p 2 . 
This is borne out by the plots of Fig. 5 in which y i2 and y 2 y are shown 
at constant fundamental voltage F, and several values of V 2 , with q> 2 
ranging ^ ip 2 ^ 2ir. This figure also establishes that y i2 and y 2x do 
not change drastically as a function of V 2 . It was also found that y l2 
and y 2i depend upon V t in an approximately linear fashion. This is 
shown in Fig. 6 where y i2 /Vi and y 2 i/Vi are plotted versus Vi for 
several values of <p 2 with V 2 constant. Thus, for moderate values of 
Vi and V 2 , we can make the approximation that y i2 and y 21 are both 
proportional to V\ and independent of s? 2 and V 2 . To demonstrate this 
analytically, let the phase of the fundamental voltage ^ 5^ 0, and con- 
sider a power series expansion of the currents i X2 = y l2 V 2 exp (j<p 2 ) and 
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Fig. 3— The input admittance, Y in i, at 3, 4, 5 and 6 GHz as it is modified by 
the presence of a second-harmonic voltage for V\ = 10 volts, V 2 = 1 volt and 
J = 340 A/cm 2 . 

*ai = 2/21^1 exp OVi). Selecting the lowest-order terms having the 
appropriate frequencies, we find that 

y l2 a T, exp (— j<pi) 



and 



(3) 



y 2l ex Vi exp OVi) 



which confirms the approximate linear dependence on V t and gives 
the appropriate form of the ipi dependence. It will be convenient later 
to approximate y l2 and y 2i by the quantities 



y l2 - K x V l exp (—#>,) = ^Vi exp [— ;(<Pi _ ML 
2/21 = -^2^! exp OVi) = k 2 V x exp [j'(vi + ML 



(4) 
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where kj = | Kj |, k 2 — \ K 2 \, <Ai = arg (ifi) and \j/ 2 = arg (K 2 ). Note 
that for pi = (only the phase <^ 2 — <Pi is important), ^ = arg (y-i 2 ) 
and 1^2 = arg {y 2l ) which is what will usually be assumed. 

The quantities y l2 and y 2x may be denned as the average of y l2 and 
y 21 over the phase <p 2 . For the 6-GHz oscillator example, the calcu- 
lated values of y l2 and y 2x as a function of frequency are shown in Figs. 
7 and 8 and the phases \p, , \p 2 and \f/ 1 + \f/ 2 are shown in Fig. 9. Obvi- 
ously these are only first-order approximations, but the complexity of 
the stability analysis requires some suitable approximation to obtain 
qualitative understanding. 

The interaction of the diode equivalent circuit of Fig. 1 with an 
external circuit can be visualized by connecting an admittance Y 2 to 
the second-harmonic port. The fundamental input admittance is then 
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Fig. 4 — The input admittance, F,„ 2 , at the second harmonic of 3, 4, 5 and 6 GHz 
for the same conditions as Fig. 3. 
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Fig. 5 — Complex plane plot of 7/12 and y 2 i for Vi = 10 volts and V 2 = 2, G, 10, 14 
volts at 6 GHz, showing the relative insensitivity of 2/12 and y 2i to changes in V 2 
and ip 2 for moderate values of V*. 



Y inl = y u — 



2/127/21 



7/22 + F 2 



(•3) 



Tuning the second harmonic by adjusting Y 2 provides the possibility 
of almost any input admittance Y inl . In particular, | Y a \ = " gives 
the short-circuit termination and Y inl = y u . Equation (5) also pre- 
dicts a pole in Y inl at the frequency for which y 22 + Y 2 = 0. This is 
not an ordinary pole as in linear circuit theory however for two reasons: 
(i) y 22 may have a negative real part because it is an active device, and 
(it) y aa is a function of V 2 so that the "pole" at y 22 + Y 2 = moves 
with changing V 2 . This means that a resonance type of behavior should 
be observed, but that the only condition where y 22 + Y 2 = is for 
Vi m 0, which is just the single-frequency oscillator condition at 2/. 
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III. STABILITY OF THE OSCILLATION STATE 

Given an oscillation state which prescribes the admittances at the 
two frequencies, there are two requirements on the circuit that must 
be met in order that this be an obtainable state of steady oscillations. 
These are the requirements of circuit realizability and oscillation-state 
stability. The realizability criterion is simply that the required circuit 
have admittances whose real parts are greater than zero. The stability 
criterion is that any perturbation away from the given state will asymp- 
totically return to the original state. 

The stability problem has been recently discussed by Kurokawa 9,10 
for the single-frequency negative-resistance oscillator. By following the 
approach used by Kurokawa and extending it to two-frequency inter- 
actions, the equations governing the stability of the harmonically 
tuned oscillator are derived in Appendix A. In this section, they are 
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Fig. 6 — Complex plane plot of K x = yn/Vj and K 2 = yn/Vi as a function of Vi 
for various values of second-harmonic phase <pi, at 4 GHz. 
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Fig. 7 — Complex plane plot of y~n and 2/21 as a function of the fundamental 
frequency for Vi = 10 volts and F 2 = 1 volt. 
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Fig. 8 — |j/i 2 | and \y tl \ versus frequency for V, = 10 volts and V 9 = 1 volt. 
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Fig. 9 — The arguments if/i, ^2 and ^i + xp* versus frequency, showing a nearly- 
linear dependence. 



applied to several special cases, and theoretical examples of their use 
with the 6-GHz germanium oscillator model of Blue are given in 
Section IV. 

In Appendix A, it is shown that the stability of an oscillation-state 
for small perturbations is determined by the solution of the system of 
equations 



whore the vector e is denned as 



+ Be = 

8a,/ V, 

8a 2 /V 2 

S(fp 2 - 2^)J 



(6) 



(7) 



and the matrix B is given by 
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As discussed in the Appendix, 5a! , 5a 2 and 5(<p 2 — 2<p x ) are the 
perturbations in the fundamental and second-harmonic voltage ampli- 
tudes and the relative phase, respectively. V t and V 2 are the unper- 
turbed values of fundamental and second-harmonic voltage amplitudes. 

The remaining quantities in the B matrix are defined as follows. The 
fundamental and second-harmonic external circuit admittances are 
Fi(w ) = G 10 + J-Sio and F 2 (2w ) = G 20 + jB 20 , respectively. The 
primes on Yi and Y 2 in equation (S) denote differentiation with respect 
to frequency at w and 2co respectively. «i and k 2 are defined in equation 

(4). 

The saturation parameters s, r and u, v are defined by equations 
(55) through (58) in the Appendix. They relate to the nonlinear satura- 
tion of the diode's conductance and susceptance at the fundamental 
and second harmonic frequencies, respectively. The significance of s 
and r is shown schematically in Fig. 10, with u and v interpreted by a 
similar diagram for the second-harmonic admittance. 

We have also introduced the angles ai and a 2 which give the slope 
on the complex plane of the circuit admittances at w and 2co 

oosa ' = vw+m • smai = vw+m (9) 



COS a, = 



(ion • ^*20 



sin a 2 = —===== (10) 



VgS + bZ' VgSTbZ 



and the angles y x and y 2 which measure the slope of the admittance 
curves y n (Vi) and y 22 (V 2 ); 

S T 

cos 7i = /-, , , i sin 7, = , ; (11) 



V ? "^T7 ' " V7+ 



U ■ V /"ION 

COS 72 = / 2 , 2 - , SU1 7 2 = / 2 2 - (12) 

Also, to and 2O are defined as in equations (48) and (49) of the Appendix 
but with the phase tp x set to zero. That is 



(13) 





0io — — <^2 — ^1 


and 






020 = <P2 — ^2- 


Note that 






010 + 020 = —^1 — ^2- 



(14) 
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Fig. 10 — Interpretation of the saturation parameters s and r for the fundamental 
admittance yu- Similar definitions hold for u and v for the second-harmonic ad- 
mittance 7/22. 



For the Ge oscillator considered here, the direct relationship between 
<p 2) 020 and 0io, as determined from equation (13), is shown in Fig. 11 
for several frequencies. 

The angles a l} 7i and 1O are shown in Fig. 12 which is a plot of the 
negative of an assumed circuit admittance —Y^a) and the diode 
single-frequency admittance 2/n(Fi) in the neighborhood of the funda- 
mental frequency. The point of intersection at o s gives the frequency 
and amplitude of the fundamental oscillation with zero second-har- 
monic voltage. As the voltage V 2 is increased by presenting an appro- 
priate value of F 2 (2w )> the frequency will shift to some new value 
o> generally accompanied by a change in voltage to V t . This shows 
that the current injected into the fundamental circuit by the ?/ 12 F 2 
expOVs) current source of Fig. 1 is just that sufficient to obtain the 
difference between the admittances — F^wo) and yn(Vi). This addi- 
tional admittance may be considered as a vector pointing from ijn(Vi) 
to — Fx(g> ), and it is the angle 1O measured clockwise about the yu(Vi) 
point that determines the orientation of this vector. Its length is given 
by 1 2/i z | ^V^i • The angle c^ gives the slope of the circuit curve at 
— Fi(a>o), and the angle 7, gives the relative change in reactive to real 
part of yu(Vi) with increasing voltage Vt at the operating point. The 
angles a 2 , y 2 and 2o may be defined in a similar manner in the second- 
harmonic admittance plane. 

The solution of equation (6) subject to a small initial perturbation 
has a decreasing amplitude with increasing time if the eigenvalues of 
the stability matrix B all have real parts greater than zero. Suitable 
tests have been devised to determine this property. 11 The general case 
is difficult to do analytically and generally difficult to interpret if done 
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numerically because of the large number of parameters of the system. 
This is done however for the 6-GHz oscillator example given in Section 
IV, and the results are compared with the simplified results of this 
section. 

In the remainder of this section, three special cases are examined 
which are severe approximations to the general case, but which yield 
interesting information. The first of these is that of a single-frequency 
oscillator, V 2 = 0. The second is the fictitious weak-coupling case 
which does not apply to the germanium diodes modeled here, but is 
included because of simplicity and for completeness. The third case is 
that of a strongly coupled small-signal approximation which gives 
qualitatively most of the features observed from the complete study 




W* 



Fig. 11 — Oscillator phase relations for the 6-GHz germanium example; ip« versus 
2o with loci of constant d l0 at 3, 4, 5, 6, 7 and 8 GHz. 



1792 



THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1970 



y„(a,) 



a =o 




-Yaw 



(-) CONDUCTANCE 

Fig. 12 — An assumed fundamental admittance plane plot showing the angles 
ai, 7i and 0io. The device admittance is yn(oi) and the negative of the circuit ad- 
mittance is — Yi(a). A similar diagram defines a 2 , 72 and 2 o in the neighborhood 
of 2 wo- 
of the eigenvalues of B, which is carried out in Section IV for the ger- 
manium diode case. 

3.1 Single-Frequency Limit 

In the very special case of V 2 = 0, only the first and third parts 
of equation (6) remain and they give the conditions 

sin(a! - Ti) > (15) 

and 



T ^ i T sin(a 2 + m ) <0. 

J 2 



(16) 



These are simply the conditions required for stability of a single- 
frequency oscillator [equation (15)] with the added condition (16) due 
to the coupling to the harmonic. If the coupling to the harmonic, 
k 2 , is zero for V 2 = 0, equation (16) does not apply. Thus, for the 
single frequency oscillator with V 2 = 0, the familiar stability relation is 
recovered. 10 
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3.2 Weak-Coupling Limit 

For an oscillator having very small k, and k 2 , the first two parts of 
equation (6) decouple. This gives 

sin(a, - 7 y) > (j = 1, 2) 

which are the single-frequency stability conditions at co and 2co for 
; = 1 and 2, respectively. The third equation then requires 

sin(a, + 1O ) + m sin(a 2 + 2O ) < (17) 

where the parameter n is defined by 



M = o 



k 2 V, | Y[ 



(18) 



2k, V 2 I Y 2 | 
We may write equation (17) as 

sm( Va + g) < (19) 

where £ is defined by the equations 

psin£ = — [sin(^, - a,) + m sin(^ 2 — a 2 )] (20) 

and 

p COS£ = — [cos(^, — a,) — M C0S(^ 2 — a 2 )]. (21) 

For a given pair of V ,, V 2 and for a fixed circuit, equation (19) thus 
gives the range of <p 2 for stable operation in the weak coupling limit. 

3.3 Small-Signal, Strong-Coupling Limit 

For very small signals the admittances y u and y 22 are independent 
of F, and V 2 so that s = r = u = v = provides another approxima- 
tion of some interest, providing that the coupling is still significant. 
In this limit, we obtain four constraints which are necessary and 
sufficient 11 to insure that the matrix B have positive eigenvalues. 
These are 

A-, = -sin (a, + 0,o) - M*m (a 2 + 0„ n ) > 0, (22) 

k 2 = -sin (a, + 0, o )-sin (a-, + 0,„) 

+ 3 cos (a, + O-cos (a 2 + 2U ) > 0, (23) 

A-., = sin (« a + 2O ) + n sin (a, + 0,„) > 0, (24) 

k, = k,k 2 - A-3 > 0, (25) 

where n is defined by equation (18), 
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The significance of this case is that for it = 1, conditions k v > 
and /c 3 > are contradictory. This implies that n = 1 is a critical 
value and is indeed unstable, whereas for n approaching zero or in- 
finity stable states of oscillation do exist. These \i « 1 and n » 1 
stable states are exclusive of each other so that, as the conditions of 
oscillation are changed, if y. passes through the value unity a discon- 
tinuity in the oscillation will occur wherein the phase, the power and 
the frequency may all jump suddenly to new values. 

To demonstrate the existence and exclusive nature of the n « 1 
and n » 1 limits, consider equations (22) through (25) Note first of 
all that if a solution is obtained for a given value of it, the solution for 
the reciprocal of that value of /x is obtained by interchanging the sub- 
scripts 1 and 2 on the angles a and 0. Thus, we need only consider the 
limit n « 1; the limit /x » 1 being obtained from symmetry. For /x « 
1, equations (22) and (24) yield [Using equation (13)] 

7T - fc + a, < tp 2 < 2tt - rP l + a, (A:, > 0) (26) 

and 

i£ 2 - a 2 < <p 2 < r + t 2 - a 2 (k 3 > 0) (27) 

respectively. For purposes of illustration we consider ai = a 2 = a. 
Then the regions defined by equations (26) and (27) may be plotted 
in the <p 2 , a plane. From equation (25), if k u /c : , and k t are > 0, k 2 > 
is automatically satisfied. Consider the constraint fc 4 > 0, which may 
be written 

— cos(a! 4- 0, o )[2 sin(a t + 6 in + a 2 4- 6 20 ) 

+ sin( ai + 6 l0 - a 2 - B 2n )] > 0. (28) 

We see that cos (a, 4- 0io) = is a critical condition, on either side 
of which the term in the brackets must also change sign. Thus, the 
lines 

<p 2 = fc - a ± tt/2 (k t = 0) (29) 

in the <p 2 , a plane are critical lines. Further, consider cos(ai + O 10 ) > 
0, then 

sinty, + ^ 2 - 2«) > -sin (2^ 2 + f , - &,)/2 (fc 4 > 0). (30) 

Equation (30) represents a curved boundary in the (p 2 , a plane and 
must be computed numerically. In Fig. 13 the regions bounded by 
equations (26), (27), (29) and (30) are plotted. The data used for this 
figure (\J/ 1 and ^ 2 ) were taken from the Ge IMP ATT example at a fre- 
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Fig. 13 — Regions of stable <pi versus a(ai = a-i = a) in the strongly coupled 
small-signal limit at 4 GHz; n « 1. 



quency of 4 GHz from Fig. 9. Figure 13 shows that, for /x « 1, there 
arc two disjoint regions. Also indicated are the values of <p 2 for which 
0io = 0, 7r/2, t, Bir/2. The angle 1O (Fig. 12) measures the relative 
location of the diode's actual input conductance with respect to the 
single-frequency large-signal negative conductance, at the fundamental 
frequency. For —ir/2 < 1O < tt/2, cos 1O is positive and the input 
conductance is less negative than it would be for zero harmonic voltage. 
For this range of 1O then, the fundamental output power is degraded 
by harmonic tuning. On the other hand, for tt/2 < IO < 3ir/2, the 
input conductance is more negative than for V 2 = 0, and the funda- 
mental output power is enhanced by the presence of harmonic tuning. 
These relationships can readily be seen by rewriting equation (1) 



Re(F inl ) = -fl,+ I 2/12 



V 2 COS 0io 



Indeed, 0, o = tt maximizes the fundamental output power for the 
particular values of V u V 2 being studied. We see that at 4 GHz, the 
maximum fundamental power point exists within a stable region for 
H « 1. It is also interesting that the minimum fundamental power 
phase (0,o = 0) is in a separate region which requires a considerably 
different circuit. 
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To obtain the similar diagram for /x » 1, the same considerations 
can be reapplied to fc x through k i} or the subscripts on <p and a can be 
interchanged. Either way, Fig. 14 shows the result. Comparison of 
Figs. 13 and 14 shows indeed the disjointed, mutually exclusive be- 
havior of the /x « 1 and y. » 1 regions of stability. Additionally, it 
shows that for a given circuit (i.e., a given a), there are two stable 
ranges of phase <p 2 (if any at all) depending on the value of m relative 
to unity. One of these encompasses the 1O = tt maximum power phase 
and the other encompasses the 1O = minimum power phase. A change 
in the bias current, which does not alter significantly the circuit vari- 
able a, may well change the relative value of n from >1 to <1 or 
vice versa, and such a change would necessitate a change of phase 
to a different branch. Thus, which branch of the stability diagram 
the oscillation state is in is determined by the history of tuning and 
bias current changes. This type of behavior would be observed experi- 
mentally as a hysteresis in frequency or power or both, which if analyzed 
would indicate that the input admittance of the diode at the funda- 
mental frequency is a nonunique function of the fundamental RF volt- 
age. The presence of this effect would be indicated if one were able 
to obtain two different values of power output for the same frequency 




Fig. 14— Regions of stable <?■> versus a(<*i = a* = a) in the strongly coupled 
small-signal limit at 4 GHz; it X> 1. 
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by changing the bias current only, without retuning the RF circuit 
in any respect. Observation at a single frequency is required in order 
to rule out the possibility of multiple-valued circuit admittances. 

In the next section, we compute the regions of stability for the 
germanium IMPATT example in full generality; that is, we use the 
complete form of the matrix B, equation (8). This must be done nu- 
merically so a limited number of cases can be examined, and the results 
are compared with the approximate forms of this section. 

IV. 6-GHZ GERMANIUM OSCILLATOR EXAMPLE 

Using Blue's approximate large-signal analysis, 1 the equivalent cir- 
cuit parameters of Fig. 1 have been calculated for a germanium diode 
of depletion layer width 4.75 microns with an assumed avalanche zone 
width of 1.5 microns. This gives a critical field E c = 1.87 X 10 5 V/cm 
for a bias current density J = 340 A/cm 2 , which agrees quite well 
with the value obtained from a more exact numerical treatment. The 
design of this model was an attempt to model the germanium diodes 
reported by Swan 2 and by Gewartowski and Morris. 12 Because the 
Read theory is slightly incorrect in its reactive effects, the frequency 
of maximum negative conductance was at about 6 GHz for the model 
but appeared to be at about 8 or 9 GHz for the actual diodes. In com- 
paring the results of this work with those of the experiments, it there- 
fore seems most useful to discuss frequency relative to f max , at which 
maximum output power is obtained. Thus, 4 GHz in this analytical 
work is roughly equivalent to 6 GHz in Swan's experiments. Table I 
lists the large-signal information obtained from Figs. 2, 8 and 9 that is 
needed for the solution of the stability constraints. This information 
was obtained for Vj = 10 volts and V 2 = 10 volts, and a dc bias current 
density J = 340 A/cm 2 . 

It is known that at resonance in a low-loss circuit where the real 
part of the admittance is constant or nearly so, the external Q can be 
written 

u dB 
WeKt 2G n da> 

where G is the real part of the admittance at co and B is the suscep- 
tance. Resonance is defined by the vanishing of J5(co ). It is useful 
here to extend this definition to define the slope parameters 

T) . = ^q dY-j. 

2G in do) ,,,.„„ 
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Table I — Diode Large-Signal Parameters At 7, = V 2 = 10 Volts 





3 GHz 


4 GHz 


5 GHz 


6 GHz 


0ii (mhos /cm 2 ) 


5.8 


12 7 


15.4 


15.6 


(722 (mhos /cm 2 ) 


15.6 


12.7 


8.5 


4.3 


—^- (mhos/cm 2 -volt) 
aV i 


0.0 


0.22 


0.24 


0.21 


—^r- (mhos/cm 2 -volt) 
av i 


1.1 


0.65 


0.30 


0.20 


— Ijp- (mhos/cm 2 -volt) 

OK 2 


0.21 


0.125 


0.065 


0.035 


— ^|- (mhos/cm 2 -volt) 
av 2 


0.20 


0.0 


0.0 


0.0 


\f/i (t radians) 


-0.2089 


-0.3056 


-0.395 


-0.477 


\f/2 (it radians) 


0.9031 


0.7742 


0.6181 


0.4798 


ki (mhos/cm 2 -volt) 


1.25 


0.86 


0.59 


0.42 


k 2 (mhos/cm 2 -volt) 


1.25 


0.70 


0.385 


0.205 



at the fundamental frequency and 



D 9 = 



G20 


dY 2 



at the second harmonic. If, at co = w and w = 2co , G' 10 and G' 2Q vanish 
respectively, then D t and D 2 reduce to the external Q's of the circuit 
at these two frequencies, particularly since the major portion of the 
diode's susceptance is considered to be part of the external circuit. 

Since, at an equilibrium point, from equations (44) and (46) of the 
Appendix 



and 



G 10 = 0i — k, V 2 cos 0j 



#20 = 02 — k 2 V x cos 2 , 



specification of the parameters D t and D 2 permits the calculation of 
I Y[ I and | Y 2 \ from the information of Table I. 
The general stability criteria for the matrix B are as follows: Let 
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B be represented 



B = 



(1 


b 


c 


(/ 


e 


/ 


.9 


h 


i 



The condition that the eigenvalues of B all be positive implies that 

fc, = a + e + i > 0, 

k 2 = ae + et + az — bd — jh — gc > 0, 

fc 3 = det B > 
and 

fc< = fc,fc a - fc 3 > 0. (31) 

These conditions must be checked numerically, and the number of 
independent variables for a general study is quite large. In the calcula- 
tions done here, the circuit variables have been restricted to a x = a 2 = 
a, with two sets of slope parameters; (i) Di = 50, D 2 = 500 and (it) 
Di = 50, D 2 = 10. The restriction on ai and a 2 is quite artificial but 
allows comparison with the approximately determined regions of Section 
III. The two sets of slope parameters D i} D 2 are an attempt to model 
(i) a high Q and (*i) a low Q second-harmonic circuit, respectively, 
and to thereby approximate the two conditions n <<C 1 and n y> 1 for 
the same set of diode data. 

The results of these calculations arc shown in Figs. 15 and 16 for 
the frequencies 3, 4, 5 and 6 GHz. These show the values of stable 
second-harmonic phase <p 2 as functions of the circuit angles, a r = a 2 = 
a. These regions repeat themselves with a periodicity of 2ir in both 
v? 2 and a. Only the principle branches are shown but it should be under- 
stood that wherever one of these regions extends across the boundaries 
chosen, it should be reflected back into the region at the opposite 
boundary- Figure 15 is for the case D x = 50, D 2 = 500, and corresponds 
to a value of n < 0.4 everywhere. Figure 16, for which D x = 50, D 2 = 
10, corresponds to values of n from near or slightly less than unity, to 
greater than 4 to 8 (the only exception is in Fig. 16a where one region 
appears having a value of n ~ 0.02). It should be noted that the value 
of n = 1 is no longer a critical value, inasmuch as stable states may 
now exist for which ft — 1. They do not appear to be large in number, 
however, and one may think of n = 1 as a transition value for which 
the area of the stable regions in the tp 2 , a plane becomes small. 
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Fig. 15 — Large-signal regions of stable ipi versus a(ai « oj = a) as obtained 
from the eigenvalues of the complete B matrix for the germanium oscillator example 
at (a) 3 GHz, (b) 4 GHz, (c) 5 GHz, and (d) 6 GHz; circuit variables Z), = 50, 
D 2 = 500; diode variables Vi = V 2 = 10 volts, J a = 340 A/cm 2 . This figure has 
li < 1 everywhere. 

Consider the 4-GHz results and compare Figs. 15b and 16b with 
Figs. 13 and 14. The locations of the stable regions in the <p 2 , a plane 
show a one-to-one correspondence but with greatly distorted shapes. 
It therefore appears that the strongly coupled small-signal approxima- 
tion used in Figs. 13 and 14, together with the n « 1 and y. » 1 cases, 
does give useful information about the general location of these stable 
regions for more realistic cases. The general properties of disjointedness 
and mutual exclusiveness are no longer strictly true (for example, 
there is some overlap of the regions centered at a = t in Figs. 15d 
and 16d). However, it is easy to see that tuning discontinuities may 
still occur, and that the circuit angles a must be considerably different 
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to obtain oscillation at 1O = ir, for example, for the two different sets 
of values of slope parameters considered. 

It is interesting that the angles a x and a 2 (and therefore, a) are 
equal to ir/2 for simple shunt resonant circuits at both co and 2cu, and 
that the stability diagrams show no cases of stable operation for this 
condition. Because of the approximations of this analysis, this cannot 
be construed to be a general conclusion, even for the diode modeled. 
It does show however, that such conditions may arise and that obtaining 
just the correct phase relations for maximum output power with a 
given circuit may be extremely difficult. 
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Fig. 16-^-Large-signal regions of stable ^ 2 versus a(an = a t = a) as obtained 
from the eigenvalues of the complete B matrix for the germanium oscillator example 
at (a) 3 GHz, (b) 4 GHz, (c) 5 GHz, and (d) 6 GHz; circuit variables Di = 50, 
D t = 10; diode variables F, = F, = 10 volts, J = 340 A/cm 2 . This figure has 
n > 1 everywhere except as noted. 
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Another observation is that the angle 6 l0 = ir for maximum funda- 
mental power output does have a stable realization in almost every 
case examined, even with the restriction a x = a 2 . 

If the points of operation along the circuit admittance curves Fi(co ), 
F 2 (2co ) are near minima of their real parts, the angles c*i and a 2 are 
restricted to lie in the range < a, < it, j = 1, 2. Such a limitation 
seems to imply different possibilities at the four frequencies calculated. 
At 3 GHz. stability is obtained in the neighborhood of d 10 = v and 
only for the D 2 = 500 case (m < 1). At 4 GHz, stability near 6 l0 = ir 
is only obtained for the D 2 = 500 (m < 1) case, but there are additional 
stable states at or near 1O = for both the D 2 = 500 (m < 1) and D 2 
= 10 (m > 1) cases. Also, at 4 GHz, Fig. 16b shows a region which 
encompasses the 8 10 = x/2 point which is a crossover between enhanced 
and degraded fundamental power. The 5-GHz cases are very similar 
to those at 4 GHz except that there are more enhanced-power stable 
states for the D 2 = 10 (n > 1) case than at the lower frequencies. 
At 6 GHz, this shift is more advanced with roughly an equal number 
of stable states in the enhanced power region for the D 2 = 10 (m > 1) 
and D 2 = 500(m < 1) cases. 

V. SUMMARY AND CONCLUSIONS 

An analysis of the stability of the tuned-harmonic mode in IMPATT 
oscillators has been presented using a simplified model of the frequency 
conversion in the avalanche diode. It has been shown that the stability 
constraints are generally quite restrictive and difficult to satisfy, par- 
ticularly for diodes showing strong harmonic interactions. The goal 
of this work has not been to present a set of design curves which insure 
stable tuned-harmonic operation, but rather to consider the difficulties 
which the stability constraints present. 

When the circuit restricts the voltage across the diode to be largely 
sinusoidal, this analysis reduces to that of the stability of a "single- 
frequency" oscillator. For nonzero fundamental and second-harmonic 
voltages Vi and V 2 , a characteristic parameter /x has been defined 
[equation (18)] which is dependent upon both diode and circuit char- 
acteristics and degree of excitation. The value of jt — 1 appears to be 
somewhat critical in that the stable regions for n > 1 and m < 1 are 
usually separate. Any tuning or bias changes which force /* to pass 
through unity are very likely to produce sudden changes in the output 
variables, i.e., power and frequency. For example, the single-frequency 
oscillator is destined to have y. ^> 1 because of the small value of V 2 . 



IMPATT DIODES 1803 

However, for equal 7, and V 2 and D 2 /D 1 ~ 10, y. < 0.4. Thus the 
single-frequency oscillator and the tuned-harmonic oscillator (high Q, 
2co circuit) are likely to operate in different regions of stability. 

The numerical treatment of the stability criteria have been restricted 
to the case where the circuit angles a t and a 2 are equal. Thus the 
results presented here cannot be considered complete. However, in the 
example studied, it was found that at an operating frequency two- 
thirds the frequency of maximum output power, the phase <p 2 for maxi- 
mum power is indeed stable and also corresponds to a realizable cir- 
cuit. It was also found that it is possible to degrade the output power, 
and therefore, harmonic interactions when improperly adjusted can 
severely lower a diode's output power from that which would exist 
with no harmonic voltage at all. 

As a necessary part of this instability analysis, a two-port model 
for the interaction was introduced and characterized for the 6-GHz 
germanium IMPATT model presented. This characterization illustrates 
the role of the second harmonic in introducing a "pseudo-pole" into 
the nonlinear admittance of the fundamental, and it clarifies the rele- 
vance of the single-frequency admittance plane characterization for 
the tuned-harmonic mode of operation. 

This analysis also has assumed that y 12 and y 2l may be described 
by equation (4). If, on the other hand, y 12 and y 2X are assumed constant, 
then this analysis becomes identical with that of two nonlinear oscilla- 
tors coupled through a linear circuit. That analysis can be carried 
through in the same manner as presented here. In such a case, the 
weakly coupled case becomes of considerable interest and has been 
treated by Schlosser. 13 

It is not necessary, of course, to introduce the two-port model of 
Fig. 1 at all, with its attendant assumptions and approximations, but 
it is possible to consider the perturbation of the oscillation-state di- 
rectly from the numerical solution of the IMPATT equations. This 
would be a more accurate method to pursue; however, it is felt that 
the approach presented in this paper provides insight that might be 
obscured in a more complicated approach. 
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APPENDIX A 

Derivation of the Stability Matrix 

In this appendix, the stability of the oscillation-state is considered 
using a linearized perturbation treatment about any general large- 
signal operating state. The result of this appendix is the derivation of 
the state-equation (6) and the stability matrix B, equation (8). 

Consider a prescribed state of oscillation satisfying the two condi- 
tions 

YM + Y**(y* » ^i . ft * ft) - (32) 

and 

F 2 (2a>o) + r«(7, ,V 2 ,<p { , <p 2 ) = 0, (33) 

where Yx(u ) and F 2 (2co ) are the circuit admittances at co and 2a> 
respectively. An approximation is made that the input admittances 
of the diode, F, nl and Y in2 , are slowly varying functions of frequency 
as compared with the circuit admittances Fi(« ) and F 2 (2w ). This is 
facilitated by considering the depletion layer capacitance, for example, 
to be a part of the external circuit. Generally speaking, equations 
(32) and (33) prescribe a functional dependence of u, the frequency 
of oscillation, upon the voltage amplitudes and phases for small varia- 
tions. For small variations in « we can approximate 



F,(co + 5,) « YM + ^ 

and 

dY, 



5, 



r 2 (2o> + « 2 ) « F 2 (2co ) + 



tlui 



d, . 



The 8 k can be determined by allowing the voltage amplitudes and 
phases to be slowly varying functions of time 

Vt(t) = a,(0 cos [oof + *>,(«)] (34) 

and 

v 2 (t) = a 2 (t) cos [2u t + <p 2 (t)]. (35) 

Differentiating with respect to time gives 
| =Re {L + 3 | i + I| 1 ] o ,exp [ ,W + „)]} (36) 
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and 

Thus, we can identify 10 



«2 cxp Li(2w f + <p 2 )] 



(37) 



and 



and therefore 



dcpi _ . 1_ dth 
1 " = dt J a, dt 



d<p 2 _ ■ 1_ ^ fl 2 
d2 ~ dt ~ 3 a 2 dt ' 



K,(« + 5,) » K.Ccoo) + 



(/a.- 



and 



F 2 (2a, + 5 2 ) » F 2 (2a) ) + ^ 






(38) 



,/<-'=£» :i!n 



are the circuit admittances related to slow variations of the amplitudes 
and phases. 

From the equivalent circuit of Fig. 1, the currents at the fundamental 
and second harmonic are 



h(t) = Re {[y u O| exp (#>,) + y l2 a 2 exp (;^ 2 )]-exp (ju t) 



and 



i 2 (t) = Re \[y 2i ai exp (#>,) + y 22 a 2 exp (jp 2 )]-exp (j2w 0l . 
which may be rewritten using the assumptions (4) as 
*'i(0 = [—Qi&i + ^a,^ cos(2y>i — <p 2 — ti)] cos(a> £ + Vi) 

+ [ — 6,o, + Kia,a 2 sin(2^, — p 2 — <£i)l sin(w £ + p,) (40) 
and 

*■(*) = [-f/ 2 a 2 + K 2 a, cos(^ 2 - 2^>j - ^ 2 )] cos(2w $ + <? a ) 

+ [-M 2 + K 2 a\ sin(v? 2 - 2^, - \f/ 2 )\ sin(2w < + ¥>»)• (41) 
Here we have introduced 

.'/u = -0i + ;'b, 
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and 

2/22 = -02 + Ah ■ 

Kirchoff's laws for the nonequilibrium case are 

i',(0 + Re { Y 1 (u> l )a l exp (jp.) exp (jw Q) = (42) 

and 

* 2 (i) + Re j F 2 (« 2 )a 2 exp (#> 2 ) exp ty2woQ) = 0, (43) 

where <a t and w 2 are the perturbed fundamental and second-harmonic 
frequencies. 

Equations (40) and (41) with (42) and (43) give the following four 
differential equations for the quantities a^t), a 2 (t), ipi(t) and <p 2 (t) 

G l - g x + G[ ^ + Bi }^= -*i<** cos «. , (44) 

-(S, + f h ) ~B[^ + Gi^= -K,a 2 sin 0, , (45) 

02 - ff 2 + G 2 ^? + Bi l^f= -K 2 a, cos 2 , (46) 

Clt d 2 0.1/ 

- (B 2 + b 2 )-B' 2 ^ + G 2 - 1 ^?= -iw, sin 2 . (47) 

Here we have denned Y x = G x + j'Bi, F 2 = (? a + i# 2 and the primes 
denote differentiation with respect to u. Also 

t = 2?, - ? 2 - ^| (48) 

and 

2 = <* 2 - 2^ - </- 2 . (49) 

Equations (44) through (47) may be rewritten so as to contain only 
a single time derivative in each 



b.'(0, - *) - Aw + w + i^rf ^ 



a, dt 
= - Kl a 2 [B[ cos 0, + Gl sin 0,|, (50) 



G?(G. " 0.) + Bf(Bi + 61) + I F{ ' 2 ' 



= -kMGI cos 0, - tffsin 0,], (51) 
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., ,2 1 da. 



B 2 (G 2 - g 2 ) - Gi(B s + b,) + | Y' 2 | 2 ^ ^ 

= -Koflj^ cos 2 + G' 2 sin 2 ], (52) 
«(<?« - g 2 ) + fl^ + 6,) + | Y 



v i2 difi 



dt 

= -k 2 Oi[^ cos 2 - ^ sin 2 ]. (53) 

Since <px is an aribitrary quantity with no physical significance, it can 
be eliminated in favor of the difference phase cp 2 — 2<px since this appears 
in both 0j and 2 . This is done by multiplying equation (51) by 2/| Y[ | 2 , 
equation (53) by 1/| Yf, | 2 and subtracting equation (51) from (53), 
giving 

dt {(P2 ~ 2vi) 



G 2 (G 2 - g 2 ) + B 2 (B 2 + y _ g^g, - gl ) + g((g, + b t ) 

I V I 2 " lv' l 2 



/G 2 cos 2 - B' 2 sin 2 \ . _ /G{ cos 0, - B[ sin 0, 
= — MM 1 y' | 2 i + 2K l a 2 \ ■ y, , 2 

(54) 

Equations (50), (52) and (54) form the set of differential equations 
for ai(t), a 2 (t) and <p 2 (t) — 2<P\{t) which will be linearized for small 
perturbations around the oscillation state. These perturbations take 
the form 

a, - I\ + 5a, , 

n> = Vt + 8a 2 , 



and 



<P2 — -<Pl = >P20 — 2<fl0 + K<P2 — 2v5l), 



where V lt V 2 , <p l0 and <p 20 are the unperturbed values of a^t), a 2 (t), 
<Pi(t) and <p 2 (t). The perturbations in the voltage amplitudes will change 
9u &i> 92, b 2 away from their values g u b 1} g 2 , b 2 which correspond to 
da-i = 8a 2 = 8(<p 2 — 2<Pi) = 0. Thus, we define the saturation parameters 
s, r, u, v which describe the linearized variation of g t around g lt etc., 
by the equations (see Fig. 10) 
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7, 5(5,0 + &,) 



^?io 



5a, 



„ = Y± *( g »° ~ </ ■■>) 

r/oo 5a 2 



and 



y = 



F 2 5(£ 20 + b 2 ) 



5a 9 



(56) 
(57) 

(58) 



where the zero subscript on the circuit variables denotes their evalua- 
tion at o> or 2a> as appropriate. 

Equations (50), (52) and (54) may now be cast in a simple matrix 
form 



dt 



+ Be = 



(59) 



where the vector e is defined as 

5a,/ 7, 1 

8th/ V t (60) 

J(<p 2 — 2<o,). 

and the matrix B is given by equation (8) of Section III. Equation (59) 
indicates that the perturbations decay with time, giving a stable state 
of oscillation, if the eigenvalues of the matrix B are all positive. 

APPENDIX B 

List of Symbols 

a t , a 2 Slowly varying amplitudes of the fundamental and second- 

harmonic voltages; equations (34) and (35). 
Stability matrix; equation (8). 

Fundamental and second-harmonic external circuit suscep- 
tances; following equation (47). 

Imaginary parts of y n and y 2 2, the susceptances of the 
single-frequency oscillator admittances; following equa- 
tion (41). 

Fundamental and second-harmonic external circuit slope 
parameters; Section IV. 

Fundamental and second-harmonic external circuit con- 
ductances; following equation (47). 



B 

Bi , B 2 



D i} D, 



6', , G s 
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f/, , r/ 2 Negative of the conductances of the single-frequency os- 

cillator admittances; following equation (41). 

K t , K 2 Complex normalized form of y l2 and y n \ equation (4). 

s, r Saturation parameters for the admittance y u ; equations 

(55) and (56). 

u, v Saturation parameters for the admittance y 22 , equations 

(57) and (58). 

F, , V 2 Fundamental and second-harmonic voltage amplitudes; 
preceding equation (1). 

F, , Y 2 Fundamental and second-harmonic external circuit admit- 
tances; Fig. 1. 

Yi n i , Y in2 Fundamental and second-harmonic IMPATT diode input 
admittances; equations (1) and (2) and Fig. 1. 

Vu , 2/22 Fundamental and second-harmonic "single-frequency" os- 
cillator admittances; Fig. 1. 

?/i2 i 2/21 Conversion transfer admittances between fundamental and 
second harmonic; Fig. 1. 

2?i2 i Z/21 Approximate form of y x2 and y 2x \ equation (4). 

a, , a 2 Fundamental and second-harmonic circuit admittance 

slope angles; Fig. 12. 

7, , 7 2 Fundamental and second-harmonic single-frequency diode 

admittance slope angles; Fig. 12. 

0! , 2 phase variables; equations (48) and (49). 

0io , #20 0i and 2 for ipi = 0, equation (13). 

k, , k 2 Magnitudes of K x and K 2 ; equation (4). 

n Stability parameter, equation (18). 

<Pi , <p 2 Fundamental and second-harmonic voltage phases; pre- 

ceding equation (1). 

\{/ i , $2 Arguments of K v and K 2 ; equation (4). 

wo Fundamental radian frequency. 
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